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Abstract. I discuss lattice QCD calculations of the properties of strongly inter- 
acting matter at finite temperature, including the determination of the transition 
temperature T^, equation of state, different static screening lengths and quarko- 
nium spectral functions. The lattice data suggest that at temperatures above 
2.0Tc many properties of the quark gluon plasma can be understood using weak 
coupling approach, although non-perturbative effects due to static magnetic fields 
are significant in some quantities. 

1 Introduction 

One expects that at sufficiently high temperatures and densities the strongly interacting matter 
undergoes a transition to a new state, where quarks and gluons are no longer confined in 
hadrons, and which is therefore often referred to as a deconfined phase or Quark Gluon Plasma 
(QGP). The main goal of heavy ion experiments is to create such form of matter and study 
its properties. We would like to know at which temperature the transition takes place and 
what is the nature of the transition as well the properties of the deconfined phase, equation of 
state, static screening lengths, transport properties etc. Lattice QCD can provide first principle 
calculation of the transition temperature, equation of state and static screening lengths (see 
ReL |l|2j for recent reviews ). Calculation of transport coefficients remains an open challenge 
for lattice QCD (see discussion in Refs. [3f?]). 

One of the most interesting question for the lattice is the question about the nature of the 
finite temperature transition and the value of the temperature Tc where it takes place. For very 
heavy quarks we have a 1st order deconfining transition. In the case of QCD with three flavors of 
quarks we expect a 1st order chiral transition for sufficiently small quark masses. In other cases 
there is no true phase transition but just a rapid crossover. Lattice simulations of 3 flavor QCD 
with improved staggered quarks (p4) using — 4 lattices indicate that the transition is flrst 
order only for very small quark masses, corresponding to pseudo-scalar meson masses of about 
60 MeV \\^ . A recent study of the transition using effective models of QCD resulted in a similar 
estimate for the boundary in the quark mass plane, where the transition is 1st order [8 . This 
makes it unlikely that for the interesting case of one heavier strange quark and two light u, d 
quarks, corresponding to 140 MeV pion, the transition is 1st order. However, calculations with 
unimproved staggered quarks suggest that the transition is 1st order for pseudo-scalar meson 
mass of about 300 MeV [11]. Thus the effect of the improvement is significant and we may 
expect that the improvement of flavor symmetry, which is broken in the staggered formulation, 
is very important. But even when using improved staggered fermions it is necessary to do 
the calculations at several lattice spacings in order to establish the continuum limit. Recently, 
extensive calculations have been done to clarify the nature of the transition in the 2+1 flavor 
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Fig. 1. The renormalized Polyakov loop Lren{T) (left) and the chiral condensate normalized to the 
zero temperature chiral condensate (right) as function of the temperature calculated at m; = 0.1ms 

[m. 



QCD for physical quark masses using Nt = 4, 6, 8 and 10 lattices. These calculations were done 
using the so-called stout improved staggered fermion formulations which improves the flavor 
symmetry of staggered fermions but not the rotational symmetry. The result of this study was 
that the transition is not a true phase transition but only a rapid crossover [12]. Even-though 
there is no true phase transition in QCD thermodynamic observables change rapidly in a small 
temperature interval and the value of the transition temperature plays an important role. The 
flavor and quark mass dependence of many thermodynamic quantities is largely determined 
by the flavor and quark mass dependence of Tc. For example, the pressure normalized by its 
ideal gas value for pure gauge theory, 2 flavor, 2-1-1 flavor and 3 flavor QCD shows almost 
universal behavior as function oi T/T^ [9j. The chiral condensate (V'V') ^-nd the expectation 
value of the Polyakov loop {L) are order parameters in the limit of vanishing and infinite quark 
masses respectively. However, also for finite values of the quark masses they show a rapid 
change in vicinity of the transition temperature. In Figure [1] I show the chiral condensate and 
the Polaykov loop as function of the temperature calculated for the pA action and light quark 
mass mi = O.lmg, with being the physical strange quark mass. The expectation value of 
the Polyakov loop vanishes in the continuum limit unless it is renormalized ( see e.g. [I3]). The 
renormalized Polyakov loop can be defined through the free energy of a static quark anti-quark 
free energy, F{r,T) normalized to the zero temperature potential at short distances |14j . 

L{T) = eM~^^), (1) 

where Foo{T) = linir^oD F{r,T). Therefore what is shown in Figure [T] is the renormalized 
Polyakov loop [TS] . In fact it is not necessary to calculate the free energy for the renormalization 
of the Polyakov loop. It is sufficient to calculate the zero temperature potential at each lattice 
spacing where finite temperature calculations have been performed [20j . This method gives 
more precise result for the renormalized Polyakov loop then the one used in Ref. [14j.. where 
the free energy has been normalized to the T = Q potential at the shortest distance. The 
difference in the value of the renormalized Polyakov loop in Figure [1] and in Ref. [20j is due to 
the different normalization of the zero temperature potential. We see that the chiral condensate 
and the renormalized Polyakov loop show rapid change at Tc suggesting that the chiral and the 
deconfinement transitions happen at the same temperature. 

To determine the value of the transition temperature and to study the interplay between 
the chiral and the deconfinement transition one usually calculates the disconnected part of the 
chiral susceptibility and the Polyakov loop susceptibility defined as 

^ = NUmf) ~ {^^n ^ - N'AiL') {Lf) (2) 
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Fig. 2. The transition temperature in units of the roTc from Ref. [7| as function of the pion mass. 



as function of the of the bare gauge coupling f3 — Here Nc^ is the spatial size of the 

lattice. The susceptibilities have a peak at some pseudo-critical coupling /3c. The chiral and 
the Polyakov loop susceptibility have been studied using lattice with temporal extent Nr = 4 
and Nr = 6 and several values of the light quark masses mi — 0.05ms, O.Itos, 0.2ms and 
0.4m<; [7J. Note that the smallest value of m; correspond to pion masses of about 140MeV. We 
find that within accuracy of the calculations pseudo-critical couplings (3c determined from the 
disconnected part of the chiral susceptibility and the Polyakov loop susceptibility coincide. This 
again shows that the chiral and the deconfinement transition happen at the same temperature. 

To determine the transition temperature we have to calculate the lattice spacing in terms of 
some physical quantity. In the past the string tension has been used to set the lattice spacing. 
A more accurate determination of the lattice spacing is provided by the so-called Sommer scale 
To defined from the static quark anti-quark potential as '^'^Jr ^ \r=ro = 1-65. Analysis of the 
quarkonium spectroscopy on the lattice lead to the value tq = 0.469(7)fm [16]. In Figure [2] I 
show the transition temperature in units of tq for different quark masses [7] and two different 
lattice spacings, corresponding to Nr — 4 and Nr — 6 lattices. Note that the value of Tc 
calculated at two different lattice spacings are clearly different. The thin error-bars in Figure 
[2] represent the error in the determination of the lattice spacing a, i.e. the error in rp/a. There 
is also an error in the determination of the pseudo-critical coupling constant f3c = 6/.g^. The 
combined error is shown in Figure [2] as a thick error-bar. For TVj = 4 calculations the error is 
dominated by the error in lattice spacing, while for A't = 6 it is dominated by the error in Pc- 
With the data on tqTc a chiral and continuum extrapolation has been attempted using the most 
simple Ansatz roTc{mTr, Nt) = roTcl^ont^'' + A(roTO,r)'' + B/Nf. From this extrapolation on gets 
the continuum value Tcrp = 0.457(7) (-1-8) (—2) for the physical pion mass m^ro ~ 0.321 7 . The 
central value was obtained using d = 1.08 expected from 0(4) scaling. To test the sensitivity 
to the chiral extrapolations d = 2 and 1 have also been used. The resulting uncertainty is 
shown as second and third error. Using the best know value of tq — 0.469(7)fm we obtain 
Tc = 192 (7) (4) MeV which is higher than the most of the previous values. It is also significantly 
higher than the chemical freezout temperature at RHIC [17]. Note that the large value of the 
transition temperature is mostly due to the large value of the string tension which is related to 
the Sommer scale as ro^/a = 1.114(4) [18]. Using the above value of ro we get y/a = 468 MeV 
which is more than 10% larger than the value ^/a = 420 MeV which was used in Ref. [19] and 
led to Tc = 173(8)MeV. Recently the transition temperature has been determined using the 
so-called stout staggered action and Nt = 4, 6, 8 and 10 [20]. The deconfinement temperature 
has been found to be 176(3) (4) MeV [20]. The central value is considerably smaller than the one 
obtained with p4 action but taking into account the errors the deviation is not very significant. 
Calculations on Nr — 8 lattices with pA action are needed to clarify this issue. The authors 
of Ref. [20] use a different definition of the chiral susceptibility which resulted in the chiral 
transition temperature oiTchirai = 151(3)(3)MeV. Using the chiral susceptibility defined above 
would result in a larger value of the transition temperature. 



4 



Will be inserted by the editor 




Lattice calculations of equation of state were started some twenty years ago. In the case of 
QCD without dynamical quarks the problem has been solved, i.e. the equation of state has been 
calculated in the continuum limit |21| . At temperatures of about ATc the deviation from the ideal 
gas value is only about 15% suggesting that quark gluon plasma at this temperate is weakly 
interacting. Perturbative expansion of the pressure, however, showed very poor convergence 
at this temperature }22j . Only through the use of new re-summed perturbative techniques it 
was possible to get agreement with the lattice data |23|24|25j . To get a reliable calculation of 
the pressure and the energy density improved action have to be used [26'27]. Very recently 
calculations with the so-called asqtad and p4 action have been done on lattices with temporal 
effects Nr — 4: and 6 [28|15j . In Figure [3] the interaction measure e — 3p as well as the pressure 
and the energy density are shown as function of the temperature for the p4 action. Calculations 
performed for — 4 and A'',- = 6 give similar results, e.g. the difference between them is at 
most 7% for the pressure. The difference in the height of the interaction measure close to Tc is 
only due to the non-perturbative features of the beta function at large lattice gauge coupling 
[15] . Therefore cutoff effects are more or less under control in these calculations. Furthermore, 
there is a good agreement between pA and aqstad calculations [28^ . We see that close to Tc the 
interaction measure is very large, which means that quark gluon plasma at this temperature is 
very far from the conformal limit, where e = Bp. At high temperature the value of the interaction 
measure is consistent with the perturbative estimate. The pressure is about 10% below the ideal 
gas value and thus it is closer to the ideal gas limit as in the previous calculations [27] . This is 
due to the fact that the previous calculations were done at bare quark masses fixed in units of 
the temperature mi/T = 0.4, while the new calculations have been performed for bare quark 
masses which correspond to fixed pion mass of about 220 MeV and physical kaon mass. This 
conclusion is also supported by the analysis of the equation of state in Ref. [29] . 



2 Spatial correlation functions 

To get further insight into properties of the quark gluon plasma one can study different spatial 
correlation functions. One of the most prominent features of the quark gluon plasma is the 
presence of chromo-electric (Debye) screening. The easiest way to study chromo-electric screen- 
ing is to calculate the singlet free energy of static quark anti-quark pair (for a recent review on 
this see Ref. [30j ) . which is expressed in term of correlation function of temporal Wilson lines 

exp(-i^i (r, T) /T) = TT:{W{r)W^ (0)) . (3) 

L = TrM^ is the Polyakov loop. Since the above correlator is not gauge invariant all calculations 
have been done in Coulomb gauge [14131132133134] . Alternatively one can insert spatial trans- 
porters between the two temporal Wilson lines which makes the correlator gauge invariant. 
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Fig. 4. The Debye mass calculated in quenched QCD [32], 2 flavor QCD [33] an in 2+1 flavor QCD. 
The lines show the leading order fit together with its uncertainty. 

Calculations of such gauge invariant correlator have also been performed and lead to results 
which are very similar to those obtained in Coulomb gauge |35|36j . The singlet free energy is 
also useful to study quarkonia binding at high temperatures [37|38|39|40|41|42j . 

In purely gluonic theory the free energy grows linearly with the separation between the heavy 
quark and anti-quark in the confined phase. In presence of dynamical quarks the free energy is 
saturated at some finite value at distances of about 1 fm due to string breaking [30|31|32 ] . Above 
the deconfinement temperature the singlet free energy is exponentially screened, at sufficiently 
large distances j32], i.e. 

Fi (r, T) = (T) - ^ ^ exp(-mz, (T)r) . (4) 

The inverse screening length or equivalently the Debye screening mass niD is proportional to the 
temperature. In leading order of perturbation theory it is mo{T) — yjl + Nf /3g{T)T. Beyond 
leading order it is sensitive to the non-perturbative dynamics of the static chromo-magnetic 
fields. The Debye screening mass has been calculated in pure gauge theory {Nf = 0) [35], in 2 
flavor QCD {Nf = 2) ^3] and very recently also in 2+1 flavor QCD [S^ and is shown in Figure 
H] at different temperatures. The temperature dependence of the lattice data have been fitted 
with the simple Ansatz motivated by the leading order result : mo{T) = + Nf /3g{T)T. 
Here g{T) is the two loop running coupling constant in MS scheme with l-tjj-g ~ 2ttT. This 
simple form can fit the data very well and we get similar values of A ~ 1.4 — 1.5 for different 
fiavor content. Thus the temperature dependence as well as the flavor dependence of the Debye 
mass is given by perturbation theory. We also see that non-perturbative effects due to static 
magnetic fields significantly effect the electric screening, resulting in about 40—50% corrections. 
However, the non-perturbative correction is the same in full QCD and pure gauge theory. Let us 
also note that in SU(2) gluodynamics the non-perturbative corrections to the Debye mass are 
approximately the same [43|44I45] . This situation can be understood in terms of dimensionally 
reduced effective theory, where the effect of hard modes with momentum p ~ ttT is integrated 
out and which contain only static electric and magnetic fields |46| . The validity of dimensional 
reduction has been tested in a wide temperature range |44I45) . 

At zero temperature the static quark anti-quark potential is determined from the Wilson 
loops: V{r) = —l/t\nW{r,t), t oo. At large separation the Wilson loop obeys the area law 
W{r,t) ~ exp(— crri) which means that the potential grows linearly with distance r. At finite 
temperature we can consider the spatial Wilson loops. They obey area law at any temperature 
Ws{x, z) ~ exp(— (Ts(T)a;2;) [47 48J. Below the transition temperature the spatial string tension is 
very close to the usual zero temperature string tension. Well above the deconfinement transition 
temperature the spatial string tension is expected to be ^/ajT) = CMg'^{T)T [iS]. This is 
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Fig. 5. The spatial string tension calculated in quenched QCD [49] and 2+1 flavor QCD [50]. The 
lines show the leading order fits. 

because in the dimensionally reduced theory it is given by CMffI and at leading order the 3- 
dimensional gauge coupling is = g'^{T)T. The spatial string tension has been calculated on 
the lattice in quenched QCD [Nf = 0) [49] and 2+1 flavor QCD [50] and the results are shown 
in Figure[5j The lattice data can be fitted very well with the simple form : yj Gg (T) — CMd^ {T)T. 
Here again g{T) is the 2-loop running coupling. For the fit we get the value of cm which agrees 
reasonably well with the result of dimensional reduction [50] . The 3-dimensional gauge coupling 
(jI has been calculated more systematically in perturbation theory and also led to a very good 
agreement with the lattice data [51I50| . 

3 Spectral functions 

Information on hadron properties at finite temperature as well as transport coefficients are 
encoded in different spectral functions. In particular the fate of different quarkonium states 
in the quark gluon plasma can studied by calculating the corresponding quarkonium spectral 
functions. On the lattice we can calculate correlation function in Euclidean time. The later is 
related to the spectral function via integral relation 



Given the data on the Euclidean meson correlator G(t, T) the meson spectral function can be 
calculated using the Maximum Entropy Method (MEM) [52] • For charmonium this was done by 
using correlators calculated on isotropic lattices [53154] as well as anisotropic lattices [55156157] 
in quenched approximation. It has been found that quarkonium correlation function in Eu- 
clidean time show only very small temperature dependence [54|57] . In other channels, namely 
the vector, scalar and axial- vector channel stronger temperature dependence was found 54 57J, 
especially in the scalar and axial-vector channels. The spectral functions in the pseudo-scalar 
and vector channels reconstructed from MEM show peak structures which may be interpreted 
as a ground state peak [55156154] . Together with the weak temperature dependence of the corre- 
lation functions this was taken as strong indication that the IS charmonia {rjc and J/'0) survive 
in the deconfined phase to temperatures as high as 1.6Tc [55|56|54j . A detailed study of the 
systematic effects show, however, that the reconstruction of the charmonium spectral function 
is not reliable at high temperatures [57]. In particular the presence of peaks corresponding to 




cosh(w(T- 1/2T)) 
sinh(w/2r) 



(5) 
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bound states cannot be reliably established. The only statement that can be made is that the 
spectral function does not show significant changes within errors of the calculations. Recently 
quarkonium spectral functions have been studied using potential models and lattice data for 
the free energy of static quark anti-quark pair [42i . These calculations show that all charmonia 
states are dissolved at temperatures smaller than 1.5Tc, but the Euclidean correlators do not 
show significant changes and are in fairly good agreement with available lattice data both for 
charnionium [54 57 and bottomonium [57 58 . This is due to the fact that even in absence 
of bound states quarkonium spectral functions show significant enhancement in the threshold 
region [42] . Therefore previous statements about quarkonia survival at high temperatures have 
to be revisited. The large enhancement of the quarkonium correlators above deconfinement in 
the scalar and axial-vector channel can be understood in terms of the zero mode contribution 
^2 59J and not due to the dissolution of the IP states as previously thought. Similar, though 
smaller in magnitude, enhancement of quarkonium correlators due to zero mode is seen also in 
the vector channel [57] . Here it is related to heavy quark transport [4141] . Although the above 
mentioned lattice studies were performed in quenched approximation we do not expect the pic- 
ture to change when dynamical quarks are included in the calculations since recent calculations 
in 2-fiavor QCD show very similar temperature dependence of charmonium correlators [60] . 

The spectral function for light mesons as well as the spectral function of the energy mo- 
mentum tensor has been calculated on the lattice in quenched approximation 61 62 63 64J. 
However, unlike in the quarkonium case the systematic errors in these calculations are not well 
understood. 



4 Conclusions 

In recent years significant progress has been made in calculating bulk thermodynamic observ- 
ables as well as spatial correlation functions on the lattice. These calculations suggest that at 
temperatures T > 2.0Tc thermodynamics can be described reasonably well using weak coupling 
approaches: re-summed perturbation theory and dimensional reduction. The temperature and 
flavor dependence of static screening length is well described by perturbation theory. However, 
the value of the screening lengths is to large extent non-perturbative and influenced or deter- 
mined by the dynamics of static magnetic flelds. Furthermore, there is no evidence for the large 
value of the gauge coupling constant at scale T. Clearly more precise lattice data and further 
perturbative calculations are needed to establish the nature of quark gluon plasma at T > 2.0Tc. 
Despite signiflcant progress our understanding of spectral functions at finite temperature is still 
quite limited. 
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